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& ' Abstract 

, We present a simple, natural #P-complete problem. Let G be a directed graph, and let k 

be a positive integer. We define q(G; k) as follows. At each vertex v, we place a fc-dimensional 
complex vector x v . We take the product, over all edges (u,v), of the inner product (x u ,x v ). 
Finally, q(G; k) is the expectation of this product, where the x v are chosen uniformly and 
independently from all vectors of norm 1 (or, alternately, from the Gaussian distribution). We 
' show that q(G; k) is proportional to G's cycle partition polynomial, and therefore that it is 

#P-complete for any k > 1. 

1 Introduction 

Let x, y € C k be £;-dimensional complex-valued vectors. We denote their inner product as 

CO 
CN 



k 

{x, y) = Y^ x*Vi ■ 

Now suppose we have a directed graph G = (V,E). Let us associate a vector x v € C k with each 
vertex v, and consider the product over all edges (u, v) of the inner products of the corresponding 
vectors: 

] ( (Xu,Xv) ■ (1) 



X 

(u,v)eE 



For instance, for the graph in Figure Q] this product is 

Y\ (Xu,X v ) = (xi,X 2 ) {x 2 ,X 3 ) (X3,xi) |(X3,X 4 }| 2 . (2) 

The expectation of this product, where each x v is chosen independently and uniformly from the set 
of vectors in C k of norm 1, is a type of moment, where each x v appears with order + d° nt . It is 
a function of the graph G and the dimension k, which we denote as follows: 

q(G;k) = Exp J j (x u ,x v ) . 
(u,v)eE 
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Figure 1: A little directed graph. We use the edge labels in the proof of Theorem [TJ 

A simple observation is that q(G; k) is zero unless G is Eulerian — that is, unless = d° nt for 
each vertex v. Since x v appears in the product djf times unconjugated and d° ut times conjugated, 
multiplying x v by e td multiplies q(G;k) by e ld ^~ d ° u '. But multiplying by a phase preserves the 
uniform measure, so the expectation is zero if 7^ d° f° r an Y v - 

So, let us suppose that G is Eulerian. In that case, what is q(G; k)? Does it have a combinatorial 
interpretation? And how difficult is it to calculate? Our main result is this: 

Theorem 1. For any k > 2, computing q(G;k), given G as input, is #P-hard under Turing 
reductions. 

If we extend #P to rational functions in the natural way, then we can replace #P-hardness in this 
theorem with ^P-completeness. 

Our proof is very simple; we show that q(G; k) is essentially identical to an existing graph 
polynomial, which is known to be #P-hard to compute. Along the way, we will meet some nice 
combinatorics, and glancingly employ the representation theory of the unitary and orthogonal 
groups. 

2 The circuit partition polynomial 

A circuit partition of G is a partition of G's edges into circuits. Let rt denote the number of circuit 
partitions containing t circuits; for instance, n is the number of Eulerian circuits. The circuit 
partition polynomial j(G;z) is the generating function 



For instance, for the graph in Figure Q] we have j(G; z) = z + z 2 . This polynomial was first studied 
by Martin |10j . with a slightly different parametrization; see also [H [3l [3^ [6l [9| 1111 112j . 
Now consider the following theorem. 

Theorem 2. For any Eulerian directed graph G = (V,E), 
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(3) 




where d v denotes d] 
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d times 



Proof. Given a vector x G C k and an integer d, the outer product of x® d = x (g) ■ ■ ■ ® a; with itself 
is a tensor of rank 2d, or equivalently a linear operator on (C fc )® rf : 

(|x> 



x 9d\( x 9d 



In terms of indices, we can write 



a.\a.2—ad 



n 



x a e x /3 e ■ 



Then n^^es ( x u-,x v ) is a contraction of the product of these tensors, where upper and lower 
indices correspond to incoming and outgoing edges respectively. For instance, for the graph in 
Figure Q] we can rewrite the product [2] as 

Y\ (x u ,x v ) = \xi) (xi\l |^ 2 ) (^2!^ \x 3 <S> x 3 ) (x 3 ® x 3 |x 4 )(x 4 |^. 
(u,v)eE 

Here we use the Einstein summation convention, where any index which appears once above and 
once below is automatically summed from 1 to k. Now, since the x v are independent for different v, 
we can compute q(G; k) by taking the expectation over each x v separately. This gives a contraction 
of the tensors 



X d = Exp 



x m\( x m 



(5) 



where d = d v , over all v. 

In order to calculate Xd, we introduce some notation. Let Sd denote the symmetric group on 
d elements. We identify a permutation n € Sd with the linear operator on (C fc )® d which permutes 
the d factors in the tensor product. That is, 



7T (X\ <g) X 2 <8> • • • ® Xd) = Z,r(l) ® X 7r ( 2 ) ® • • • ® ^(rf) , 



or, using indices, 



7T 



/9i/3 2 -/3 d 



where (5* is the Kronecker delta operator, <5*- = 1 if i = j and if i 7^ j. Diagrammatically, ir 
is a gadget with (i incoming edges and d outgoing edges, wired to each other according to the 
permutation n. 

We have the following lemma: 

Lemma 1. With X d defined as in ([5]), if x is uniform in the set of vectors in C k of norm 1, then 



(fc-1)! v 
Xd= (k + d-l)\ 



(6) 
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Proof. First, Xd is a member of the commutant of the group U(k) of k x k unitary matrices, since 
these preserve the uniform measure. That is Xd commutes with U® d for any U G U(k). By Schur 
duality, the commutant is a quotient of the group algebra CfS^]; namely, the image of C[Sd] under 
the identification above. Thus Xd is a superposition of permutations, Yln€S d a "^ 7T - 

We also have Xdir = nXd = Xd for any ir. Thus Xd is proportional to the uniform superposition 
on Sd, or equivalently the projection operator II sym = (l/d\) X^ 71 " onto the totally symmetric 
subspace V sym of (C fc )® d . Since tr Xd = Exp^ \x\ = 1 while trll sym = dim V sym , we have Xd = 
(1/ dimV; ym )n sym . 

Finally, dimV^ ym is the number of ways to label the d factors of the tensor product with 
basis vectors {ei, . . . , e^} in nondecreasing order — or, for aficionados, the number of semistandard 
tableaux with one row of length d and content ranging from 1 to k. This gives dimT4 ym = 

To illustrate some ideas that will recur in the next section, we give an alternate proof. First, 
note that tr-zr is the number of ways to label each of 7r's cycles with a basis vector ranging from 1 
to k, or k°^ where c(tt) denotes the number of cycles (including fixed points). Thus 

tr E * = E kC(7T) ■ ( 7 ) 

To compute this generating function, we use the fact that each permutation ir G Sd appears once in 
the following product, where 1 denotes the identity permutation, and Tj,- denotes the transposition 
of the ith and jth object: 

7T = 1(1 + T lj2 )(l + 7-1,3 + T 2 , 3 ) ■■■(! + T M + T 2jd + • • • + T d _ M ) . (8) 

This product works by describing a permutation lit of t objects inductively as a permutation irt-i 
of the first t — 1 objects, composed either with the identity or with a transposition swapping the 
tth object with one of the previous t — 1. If we apply the identity, then the tth object is a fixed 
point, and c(irt) = c(ir t -i) + 1, gaining a factor of k in ([7]); but if we apply a transposition, then 
c(-Kt) = c(7r t _i). Thus (JEJ) becomes 

£ = k{k + l)(k + 2) ■ ■ ■ (k + d - 1) = ik + d _ ~ f . 

Comparing traces again gives ([6]). □ 

All that remains is to interpret the operators Xd v , and their contraction, diagrammatically. 
Lemma [1] tells us that, for each vertex v of G, taking the expectation over x v converts it to a sum 
over all d v \ ways to wire the incoming edges to the outgoing edges. But doing this at each vertex 
gives us a sum over all cycle partitions of G. Contracting these tensors gives the number of ways to 
label each cycle in a each partition with a basis vector ranging from 1 to k, so each cycle contributes 
a factor of k. Along with the scaling factor in ([6]), this completes the proof. □ 

Next we show that the cycle partition polynomial is ^P-hard. To our knowledge, the following 
theorem first appeared in [7]; we prove it here for completeness. 

Theorem 3. For any fixed z > 1, computing j(G; z) from G is #P-hard under Turing reductions. 
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Figure 2: A planar graph G (black) and its oriented medial graph G m (gray). 

Proof. Recall that the Tutte polynomial of an undirected graph G = (V, E) can be written as a 
sum over all subsets S of E, 

T(G; x,y) = J2( x ~ l) c(S) ~ c(G) {y - l) c ^+\ s \~ n . (9) 

SCE 

Here c(G) denotes the number of connected components in G. Similarly, c(S) denotes the number 
of connected components in the spanning subgraph (V, S), including isolated vertices. When x = y, 
we have 

T(G; x,x) = Y,(*- l) c ^+^- c ^ , (10) 

SCE 

where £(S) = c(s) + \S\ — n is the total excess of the components of S, i.e., the number of edges 

that would have to be removed to make each one a tree. 

If G is planar, then we can define a directed medial graph G m as in Figure [2j Each vertex of 

G m corresponds to an edge of G, edges of G m correspond to shared vertices in G, and we orient 

the edges of G m so that they go counterclockwise around the faces of G. Each vertex of G m has 
d in _ d out = 2j go Gm ig Eulerian _ 

The following identity is due to Martin [ID]; see also [II], or [2] for a review. 

j(G m ;z) = z c ^T(G;z + l,z + l). (11) 

We prove this using a one-to-one correspondence between subsets S C E and circuit partitions of 
G m . Let v be a vertex of G m , corresponding to an edge e of G. Then the circuit partition connects 
each of v's incoming edge to the outgoing edge on the same side of e if e G 5, and crosses over 
to the other side if e ^ S. It is easy to prove by induction that the number of circuits is then 
c{S)+£(S), in which case {ID]) yields (fTT|) . 

The theorem then follows from the fact, proven by Vertigan [13], that the Tutte polynomial for 
planar graphs is #P-hard under Turing reductions, except on the hyperbolas (x — l)(y — 1) € {1, 2} 
or when (x,y) £ {(1, 1), (— 1, — 1), (u,u*), (ui*, ui)} where oj = e 27 ™/ 3 . Thus computing j(G; z) for 
any z > 1 is #P-hard, even in the special case where G is planar and where every vertex has 
d m = d ont = 2. □ 
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3 Real- valued vectors 



We can also consider the case where the x v are real- valued, and are chosen uniformly from the set 
of vectors in R fc of norm 1. In this case, the inner product (x u ,x v ) becomes symmetric, so the 
graph G becomes undirected. We might then expect q(G; k) to be related to the circuit partition 
polynomial for undirected circuits, and indeed this is the case. 

We again wish to compute the tensor = Exp x \x® d ) First, let denote the set of 

perfect matchings of 2d objects; note that 

\M d \ = (2d - 1)!! = (2d - l)(2d - 3) ■ ■ ■ 5 ■ 3 • 1 = ^ . 

We can identify each matching // G with a linear operator on where the first d objects 

correspond to upper indices, and the last d correspond to lower indices. However, in addition to 
permutations that wire upper indices to lower ones with a bipartite matching, we now also have 
"cups" and "caps" that wire two upper indices, or two lower indices, to each other. For instance, if 
d = 2 then includes three operators, corresponding to the three perfect matchings of 4 objects: 

Wl=\ [ 6 TJt=Yy and ^ ,oa *AA=X- (12) 

The first two of these operators correspond to the identity permutation and the transposition t\ % 
respectively, as in the previous section. The third one is a cupcap; it is the outer product of the 
vector Yli= l e i ® e i with itself, where &i denotes the ith basis vector in R fe . We denote it 7^2, and 
more generally 7^ = 8 ai,a i Sp^pj. 

Now, in the real- valued case, Lemma Q] becomes the following: 

Lemma 2. If x is uniform in the set of vectors in M. k of norm 1, then 

1 ^ _ (fc - 2)!! ^ 

d ~ k(k + 2)(k + 4) • • • (k + 2d - 2) ^ ^ ~ (k + 2d - 2)!! ^ M ' ( } 

where nil = n(n — 2)(n — 4) • • • 6 ■ 4 ■ 2 if n is even, and n(n — 2)(n — 4) • • • 5 • 3 • 1 if n is odd. 

Proof. Analogous to the complex case, is a member of the commutant of the group 0(k) of k x k 
orthogonal matrices, since these preserve the uniform measure. That is, X d commutes with 0® d 
for any O G 0(k). The commutant of 0(k) is the Brauer algebra; namely, the algebra consisting 
of linear combinations of the operators fj, G M^. Thus X d is of the form YlaeM d a^fi. 

In addition to being fixed under permutations as in the complex case, Xd is also fixed under 
partial transposes, which switch some upper indices with some lower ones. Thus Xd is proportional 
to the uniform superposition Y^u&M ^° the cons tant of proportionality, we again compare 
traces. 

As in the case of permutations, the trace of an operator fi G is k c ^\ where c([i) is the 
number of loops in the diagram resulting from joining the upper indices to the lower ones. For 
instance, for the operators in (|12p . we have tr 1 = k 2 , tr ti 2 = k, and tr 71 2 = k. Thus we wish to 
calculate 
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We can write X^eAf d 

product, analogous to ([8]): 
^ vr = 1(1 + ri j2 + 71,2) (1 + n,3 + 7i,3 + t 2) 3 + 72,3) • • • (1 + T 1>d + 7i jd H h r d _i )d + Jd-i,d) ■ 

TreS d 

This product describes a matching //t of 2t objects inductively as a matching fa—l of the first 
2(i — 1) objects, composed either with the identity, or with a transposition or cupcap connecting 
the tth upper object with the ith lower one and the tth lower object with the ith upper one, or 
vice versa. If we apply the identity, then the ith upper object is matched to the tth lower one, and 
c(fJLt) = c(/i£_i) + 1, gaining a factor of k in (TBI) ; but if we apply a transposition or cupcap, then 
c(fJLt) = c(nt-i)- Thus (fT4"j) becomes 

£ fc*0 = fc(fc + 2)(fc + 4) • • • (k + 2d - 2) = {k + f U . 

irES d 

We again have irXj = Exp x \x\ 2d = 1, and comparing traces gives (fTHj) . □ 

As before, g(G; fc) is a contraction of the tensors Xj. However, now G is undirected, with no 
distinction between incoming and outgoing edges, so at each vertex of degree d v the appropriate 
tensor is Xd v / 2 - Applying Lemma [2] to each v sums over all the ways to match v's edges with each 
other, and hence sums over all possible partitions of G's edges into undirected cycles. The trace of 
the resulting diagram is again the number of ways to label each cycle with a basis vector. So, if 
define a polynomial jundirected^; z) as Ylt^i r tz l ■, where r t is the number of partitions with t cycles, 
then Theorem [2] becomes 

Theorem 4. For any undirected graph G = (V, E) where every vertex has even degree, if we define 
q(G; k) by selecting the x v independently and uniformly from the set of vectors in M fc with norm 
1, then 

q(G; k)=(j[ {k ^~ v m 2)l ^ iundirectedtG; k) . (15) 

To our knowledge, the computational complexity of j U ndirected(G; z) is open, although it seems likely 
that it is also #P-hard. 



4 The Gaussian distribution 



Our results above assume that each x v is chosen uniformly from the set of vectors in C fc or R fe of 
norm 1. Another natural measure would be to choose each component of x v independently from 
the Gaussian distribution with variance l/k, so that Expfja^l 2 ] = 1. 
For vectors in C k , the probability density of the norm |x| is then 



p(\x\) 



2k k+1 
k\ 



,2k-l e -k\x\ 2 



(16) 



Compared to the case where \x v \ = 1, each x v contributes scaling factor of |x,_,| 2rf to the product ([T]). 
The even moments of (1161) are 



Exp 



1 2d 



(rf+fc-1)! 
k d (k- 1)! 



7 



so in the Gaussian distribution @ becomes 



(17) 



where m denotes the number of edges. 

We could also have derived this directly from the Gaussian analog of Lemma [TJ If x is chosen 
according to the Gaussian distribution on C k , and we again let X d denote Ex Pa , \x® d ) (x® d \, then 



(18) 



Tves d 



Similarly, in the real-valued case, if we choose each component of i 6 1* from the Gaussian 
distribution on R with variance 1/k, then (]15p becomes 



l(G;k) — ■— iundirected ip't k) 



(19) 



since 



Xd= yt Yl 



(20) 



fieM d 



Both (jTSjl and ([20]) are consequences of Wick's Theorem [8j H5], that if xi, . . . ,X2t obey a multi- 
variate Gaussian distribution with mean zero, then 



Exp 



' it 

n 

i=l 



5^ n ex pn^-]- 

neM t (i,j)<=n 
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